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Abstract. Since their conception in 1975, Genetic Algorithms have been an 
extremely popular approach to find exact or approximate solutions to opti- 
mization and search problems. Over the last years there has been an enhanced 
interest in the field with related techniques, such as grammatical evolution, 
being developed. Unfortunately, work on developing genetic optimizations for 
low-end embedded architectures hasn't embraced the same enthusiasm. This 
short paper tackles that situation by demonstrating how genetic algorithms 
can be implemented in Arduino Duemilanove, a 16 MHz open-source micro- 
controller, with limited computation power and storage resources. As part of 
this short paper, the libraries used in this implementation are released into the 
public domain under a GPL license. 
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1. Introduction 

In the 1950s, while studying industrial organizations, Herbert Simon pointed out 
that when confronted with real-world problems, human beings have a tremendous 
difficulty with maximizing the outcomes. His rationale is that we lack the cognitive 
abilities required to evaluate all outcomes with sufficient precision, with our weak 
and unreliable memories. He stated that when present with a particular problem, 
we are forced to make decisions not by maximization but by satisficing [8 . This sat- 
isficing criterion consists of setting a particular aspiration level which, if achieved, 
they will be happy enough with, and if they don't, try to change either their aspira- 
tion level or their decision. A Genetic Algorithm (GA), is a search space algorithm, 
which follows this satisficing criterium. In most of the optimization problems, often 
times it is more important to obtain a solution that meets our requirements than 
to obtain an optimal solution. In other words, when confronted with a problem, we 
are looking for some good , satisficing level of performance quickly with reasonable 
resources. 



1.1. Genetic Algorithms. A genetic algorithm is a technique used in computing 
to find exact or approximate solutions to optimization and search problems. Ge- 
netic algorithms are a particular class of evolutionary algorithms (EA) which uses 
techniques inspired by evolutionary biology such as inheritance, mutation, selec- 
tion, and crossover. Genetic algorithms are implemented in a computer simulation 
in which a population of many individuals, each containing a potential solution to 
the proposed problem. Traditionally, these solutions are represented in binary as 
strings of Os and Is, appropriately named chromosomes. The GA algorithm starts 
with a population of randomly generated individuals, and in each generation, the 
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fitness of every individual in the population is evaluated, with multiple individuals 
stochastically selected from the current population (based on their fitness), and 
modified (recombined and possibly randomly mutated) to form a new population. 
The new population is then used in the next iteration of the algorithm. The algo- 
rithm terminates when the maximum number of iterations have been processed or 
when a satisfactory fitness level has been reached for the population. Clearly, if the 
algorithm has terminated due to a maximum number of generations, a satisfactory 
solution may not have been reached. An introductory reference to GAs can be 
found in [ri^TT. 

Genetic algorithms are suited to a variety of search problems, provided one 
he can encode solutions of a given problem into a binary string (chromosomes) 
and evaluate the quality of that particular solution. The lasting appeal of GAs 
is brought from its simplicity and ability to quickly discover solutions for difficult 
high-dimensional problems. GAs have shown to be useful when: 

• Search space is large, complex or poorly understood. 

• Domain knowledge is scarce or expert knowledge is difficult to encode to 
narrow the search space. 

• No mathematical analysis is available. 

• Traditional search methods fail. 

Genetic Algorithms have been used extensively in many diverse areas on different 
optimization problems, such as the study of social systems and cooperations [3], 
numerical optimizations .2! , VLSI cell area optimization T , quantum circuit de- 
sign [12], synthesis of QCA circuits [10], index fund portfolio selection job sched- 
uling [6] and antenna design [4] . In this paper we release an implementation of GAs 
in the Arduino Duemilanove (ATmega328 - rev 2009b), a low-end micro-controller. 

1.2. Arduino Duemilanove. Arduino Duemilanove is an open-source electronics 
prototyping platform consisting of an 8-bit Atmel AVR micro-controller clocked at 
16 MHz. It is powered at 5V and contains 32 KB of fiash memory contents of which 
2 KB used by boot-loader. These specifications should be enough to demonstrate 
the limitations of this particular micro-controller. We decided to implement a 
genetic algorithm on an arduino micro-controller for the following reasons: 

• Inexpensive and easily obtainable: on February 2020 the ATmega328 unit 
price was $4.30. 

• Open source: the IDE software and the hardware reference files have both 
been released under the GPL open-source license. 

• Development Ease: the core implementation language closely resembles 
C, and its graphical IDE is freely available on multiple popular operating 
systems. 

• Established user base: the arduino platform has been used is many intro- 
ductory embedded systems university-level classes. 

• Low power requirements: On standby mode the arduino micro-controller 
requires a current of 0.2 mA. If an average arduino application requires 
around 100mA, that would entail the power consumption to be on the 
proximity of 0.5 watt. 

This GA implementation will also serve as a proof of concept in which we demon- 
strate that even extremely limited hardware can still be used in solving practical 
real world problems. 
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2. IMPLEMENTATION 

Section|4]contains the source code for the Arduino GA hbrary, GA.cpp and GA.h, 
and a sample optimization problem, main.cpp. To run this code, one must copy 
the main.cpp source into the active Arduino project, and following the instructions 
on the Arduino support documentation include GA.cpp and GA.h as a single 
library. By running main, cpp and selecting the Serial Monitor option, the outcome 
of the sample optimization problem will be displayed. 

2.1. Sample optimization problem. The main idea behind a GA is that we 
have a particular problem we want to know a possible solution for it. We do not 
know the optimal solution a priori, but we can tell how good a particular solution 
is. A GA run starts with a certain number of individuals with all its chromosomes 
set to a random binary state. These random bits can be encoded into a solution 
for the given problem. Our sample optimization problem consists on maximizing 
the function f{x) — x on the integer interval [0,31]. In other words, we want our 
best individual to have as many bits set to logic one as possible. Each time a new 
iteration is processed, each individual will have a modified set of chromosomes, 
which will be evaluated by our fitness function. This fitness function evaluator is 
coded inside main. cpp. This fitness function can be adjusted to solve any particular 
problem. In fact, the main difficulty with GAs is the encoding of a solution in with 
a binary state. Example of binary encodings for various problems can be found 

mm 

2.2. Limitations and Discussion. This GA implementation contains just a sim- 
ple bare-bones, single crossover point as described in [l^. Due to hardware con- 
straints, these algorithms have some limitations. 

• The population size is limited to 100 individuals. This is user specified in 
main. cpp. 

• Each individual has a hardwired 32 bit chromosomes. This means that the 
solution encoding is constrained by those number of bits. 

• The value of each fitness evaluation variable, which states how good a 
particular solution is, must be between and 100. 

3. Conclusion 

Our unbridled technological lust is constantly pressuring the development of 
newer, faster and smaller micro-controllers. In our efforts to develop algorithms to 
high-end, and expensive micro-controllers, lower-end devices are being neglected. 
This particular work described an implementation of a genetic algorithms cus- 
tomized for the Arduino micro-controller. In here we demonstrated that the Ar- 
duino micro-controller embedded system can successful be used to solve search 
space problems at a low cost. Future work will consist of investigating practical 
engineering uses where GA implementations on Arduino micro-controllers can have 
a distinct impact, by taking advantage of its open-source nature, low cost and low 
power requirements. This may involve parallelizing this implementation across mul- 
tiple Arduino micro-controllers. This approach will involve storing information in a 
single arduino board, and using it to dynamically program connected blank micro- 
controllers. Further work also needs to be performed as to assert the practicality of 
using Arduino micro-controllers to solve complex engineering problems regarding 
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different metrics, such as computation time, power consumption and cost per clock 
cycle. 

4. Source code 

The source code for this GA arduino library has now been released into the 
public domain under a GPL license and can be directly copied from this paper, or 
downloaded from http://www.nunoalves.com/source/GA-revA.zip. 

4.1. main.cpp. 

//main.cpp : Sample arduino entry point tliat uses tlie arduinoGA library. 
//Copyright (C) 2010 Nuno Alves 

#include <GA.h> 

void setup () 
{ 

Serial. beginOeOO) ; 
randomSeed(analogReadCO) ) ; 

> 

int populatioiisize=99; //select size in the range [0..99] 
int nuingenerations=100; //number of generations cannot exceed 65k 
//the unit of bit mutation is per thousand. For example, if we set the 
//bitmutation to 10, approximately every 10 out of Ik bits will be mutated 

int bitmutation=l ; 

//initializing the arduinoGA library class. 

//this will initialize a set of random genes. 

GA ga(populationsize,numgenerations, bitmutation) ; 

//we start in the generation 
int generation=0; 

void loopO 
{ 

if (generation==0) { Serial. printlnC") ; Serial. printlnC") ; } 

while (generation<numgenerations) 
{ 

//—=================================================== 

//Fitness Evaluation 

//====================================================================== 

//This is the part where we "how good" a particular combination of 
//genes is. In this simple case, we just want to maximize the number of I's 
//in the genes. The I's can be for example, the number of lit LED ligths. 
//NOTE: due to memory constraints, the fitness value CANNOT(!!!) be greater 
//than 100. 

//begin: Fitness Evaluation 

for (int i=0 ; i < populationsize ; i++) 

{ 

unsigned int tO_f itness=0; 

unsigned int tO_a_population=ga.read_tO_a_population(i) ; 
unsigned int tO_b_population=ga.read_tO_b_population(i) ; 

///we are evaluating the fitness on 32 bits... dividing by 
// ga . tO_a_populat ion [i] and ga . tO_b_populat ion [i] 
for (int j=0 ; j < 16 ; j++) 
{ 

if (bitRead(tO_a_population, j)==l) 
{ 

tO_f itness=tO_f itness+1 ; 

> 



if (bitRead(tO_b_population, j)==l) 
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{ 

tO_f itness=tO_f itness+l ; 

} 

} 

ga.write_tO_f itiiess(i,tO_f itness) ; 

> 

//end: Fitness Evaluation 

//the next steps perform the genetic optimization process 
//behind the scenes. It mates the best individuals and create the 
//next generation with a new set of genes. 

ga . process_generat ion (generation) ; 

//change the reportStatistics value from to 3 for more detailed printout 
ga. report Statistics (generation, 0) ; 
ga.prepare_next_generation() ; 
generation=generation+l ; 

3" //while (generation<numgenerations) 

}//end of void loopO 



.2. GA.h. 



//GA.h : header file for the arduinoGA library. 
//Copyright (C) 2010 Nuno Alves 

#ifndef GA_h 
#define GA_h 



#include "WProgram.h" 
#define GPOPSIZE 100 



class GA 
{ 

public : 

GA(int _popsize, int _ngenerations , int _bitmutation) ; 
void process_generationCint generationid) ; 

void reportStatistics (int generationid, int verbosityLevel) ; 
void prepare_next_generation() ; 

void write_tO_f itness(int idx, unsigned int value) { tO_f itness [idx] =value ; } 
unsigned int read_tO_a_population(int idx) { return(tO_a_population[idx] ) ; } 
unsigned int read_tO_b_population(int idx) { return(tO_b_population[idx] ) ; }■ 



private : 

unsigned int tO_f itness [GPOPSIZE] ; 

unsigned int top_f itness_val ; 

unsigned int sum_f itness_val ; 

unsigned int tl_a_population [GPOPSIZE] ; 

unsigned int tl_b_population [GPOPSIZE] ; 

unsigned int bestcandidate_a; 

unsigned int bestcandidate_b; 

unsigned int tO_a_population [GPOPSIZE] ; 

unsigned int tO_b_population [GPOPSIZE] ; 

unsigned int next _gen_expected_count [GPOPSIZE] ; 



void randomize_tO_population() ; 
void evaluate_tO_f itnessO ; 
void expected_count_tl () ; 
void populate_tl () ; 
void mate_tl(); 



int bitmutation; //should be within the range [0..1000] 

int popsize; 

int ngenerations; 



#endif 
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4.3. GA.cpp. 

//GA.cpp : arduinoGA library. This file contains the function implementations 
//Copyright (C) 2010 Nunc Alves 

#include "WProgram.h" 
#include "GA.h" 

GA::GA(int _popsize, int _ngenerations, int _bitmutation) 

{ 

popsize = _popsize; 

ngenerations = _ngenerations ; 
bitmutation = _bitmutation; 



> 



//initialize the entire population 
randomize_tO_population() ; 



//verbosityLevel=0 ,1,2 

void GA: :reportStatistics(int generationid, int verbosityLevel) 
{ 

if (verbosityLevel>=0) 
{ 

Serial. print ("generation=") ; 

Serial. print (generationid) ; 

Serial. print (" , ") ; 

Serial .print ("top fitness="); 

Serial .print (top_f itness_val) ; 

Serial. print (" , ") ; 

Serial. print ("sum fitness="); 

Serial. println(suin_fitness_val) ; 



if (verbosityLevel>=l) 

{ 

Serial. printC'top candidate (BIN) : ") ; 
Serial .print (best candidate_a, BIN) ; 
Serial . print ( " " ) ; 

Serial . print In (bestcandidate_b , BIN) ; 



} 



Serial. printC'top candidate (DEC) : ") ; 
Serial .print (bestcandidate_a,DEC) ; 
Serial . print ( " " ) ; 

Serial . print In (bestcandidate_b , DEC) ; 



if (verbosityLevel>=2) 
i 

Serial. printlnC'printing each individual+f itness") ; 
for (int i=0 ; i < popsize ; i++) 

{ 

Serial . print (tO_a_population [i] ,BIN) ; 
Serial . print ( " " ) ; 

Serial . print CtO_b_population [i] ,BIN) ; 
Serial . print C " ( ") ; 
Serial .print (tO_f itness [i] ) ; 
Serial. printlnC) ") ; 

> 

Serial. printlnC'printing next gen expected count"); 
for (int i=0 ; i < popsize ; i++) 

{ 

Serial . print C"individual=") ; 
Serial . print (i) ; 
Serial. pr int (" , ") ; 

Serial .println (next _gen_expected_count [i] ) ; 

} 



Serial. println("printing each of next gen individuals"); 
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for (int i=0 ; i < popsize ; i++) 
i 

Serial. print(tl_a_populatioii[i] ,BIN) ; 
Serial . print ( " " ) ; 

Serial. printlii(tl_b_populatioii[i] ,BIN) ; 

> 

> 

> 

//processes all the steps for current generation 
void GA: :process_generation(int generationid) 

{ 

//calculates the average and find the best individual 

evaluate_tO_f itnessC) ; 

//whats the probability of a particulair individual 
//show up on tl? 
expected_count_tl () ; 

//populate the next generation array with the best performing individuals 

populate_tl() ; 

//print some pre-mate statistics 
//report_pre_mate_statistics (generationid) ; 

//mate every two elements at a random crossover point 
mate_tl() ; 

> 

//randomize all genes of each individual in the tO_population 

void GA: :randomize_tO_population() 

{ 

for (int i=0 ; i < popsize ; i++) 
{ 

unsigned int randNumberA = randomO; 
tO_a_population[i] =randNumber A ; 

} 

for (int i=0 ; i < popsize ; i++) 

■c 

unsigned int randNumberA = randomO; 
tO_b_population [i] =randNumberA; 

} 

} 

//store just the current generation fitness to save memory 

void GA: :evaluate_tO_f itnessO 

{ 

top_f itness_val=0 ; 
sum_f itness_val=0 ; 

for (int i=0 ; i < popsize ; i++) 
■C 

sum_f itness_val = tO_f itness [i] +sum_f itness_val ; 
if (tO_f itness [i] >top_f itness_val) 

{ 

top_fitness_val=tO_f itness [i] ; 

//begin: store the best candidate for viewing 

bestCcindidate_a=tO_a_population [i] ; 

bestCcindidate_b=tO_b_population [i] ; 

//end: store the best candidate for viewing 

} 

} 

} 

void GA : : expected_count_tl() 

{ 

//sort all elements in ascending order 

//based on the fitness stored in tO_f itness [i] .. . 



//begin: bubble sort the array 
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for Cint x=0; x < popsize-1; x++) 
{ 

for (int y=0; y < popsize-x-1; y++) 
{ 

if (tO_fitiiess[y] > tO_f itness [y+1] ) 
{ 

unsigned int t=tO_f itness [y] ; 
unsigned int A=tO_a_population [y] ; 
unsigned int B=tO_b_population[y] ; 

tO_fitness[y]=tO_f itness [y+1] ; 
tO_f itness [y+1] =t ; 

tO_a_population [y] =tO_a_population[y+l] ; 
tO_a_population [y+1] =A; 

tO_b_population[y] =tO_b_population [y+1] ; 
tO_b_population [y+1] =B ; 

> 

> 

> 

//end: bubble sort the array 

//since we limited each fitness element not to surpass 100 
//by letting 100 population size we can have a final sum of 10k. 
//however the unsigned integer can hold at least the number 65,535. 
//this means we could theoretically have a population size of 650. 

//next_gen_expected_count [] will have a 
//cdf on the next _gen_expected_ count [] array 
unsigned int current _value=0; 
for (int i=0 ; i < popsize ; i++) 

■c 

current_value=current_value+tO_f itness [i] ; 
next_gen_expected_count [i] =current_value ; 

y 

} 

void GA: :populate_tl() 

{ 

for Cint i=0 ; i < popsize ; i++) 

■c 

unsigned int randomValue=random(0,next_gen_expected_count [popsize-1] ) ; 
int select_index=0; 

for (int j=l ; j < popsize ; j++) 
{ 

if ( (randomValue>=next_gen_expected_count [j-1] ) ft& (randomValue<next_gen_expected_count [j] ) ) 

{ 

select_index=j ; 
j=popsize; //exit the loop 

> 

> 

tl_a_population[i]= tO_a_population[select_index] ; 
tl_b_population[i]= tO_b_population[select_index] ; 

> 

> 

//switch bits between two individuals at a particular cross point 
void GA: :mate_tl() 

{ 

for (int i=0 ; i < popsize ; i=i+2) 

■c 

int crossoverSite=random(0,31) ; 
int crossover_siteA; 
int crossover_siteB; 



//begin: mutate genes a bit. For simplicity we just randomize an entire 
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//set of genes . . . 

unsigned int randomMutationA=random(l , 1000) ; 
unsigned int randoinMutationB=random(l , 1000) ; 

if (randomMutationA<=bitiiiutation)-C tl_b .population [i] = randomO; } 
if (randomMutationB<=bitiiiutation)-C tl_b .population [i] = randomO; }■ 
//end: mutate genes a bit 

if (crossoverSite<16) 

■c 

crossover_siteA=crossoverSite; //keep values on array A till crossoversite 
crossover_siteB=0; //cross all values on array B 

> 

else 
{ 

crossover_siteA=16 ; //dont mate any values of array A 
crossover_siteB=32-crossoverSite; 

> 



//mate array A 

for (int j=crossover_siteA ; j < 16 ; j++) 

{ 

bool siteA=bitRead(tl_a_population[i] , j) ; 
bool siteB=bitRead(tl_a_population 
bitWrite Ct l_a_population [i] ,j ,siteB) ; 
bitWrite Ctl_a_population [i+1] , j ,siteA) ; 

} 

//mate array B 

for (int j=crossover_siteB ; j < 16 ; j++) 

{ 

bool siteA=bitRead(tl_b .population [i] , j) ; 
bool siteB=bitRead(tl_b_population [i+1] , j) ; 

bitWrite Ctl_b_population [i] , j ,siteB) ; 
bitWrite (tl_b_population [i+1] , j ,siteA) ; 

} 

} //for (int i=0 ; i < popsize ; i=i+2) 

} 

//this function just copies each element of tl into tO 

void GA: :prepare_next_generation() 

i 

for (int i=0 ; i < popsize ; i++) 

{ 

tO_a_populatioii [i] =t l_a_populatioii [i] ; 
tO_b_populatioii[i]=tl_b_population[i] ; 

} 

} 
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